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This paper is a development of the note [1]. The purpose of the 
paper is to deduce a number of equations which describe excitations 
propagating along n- dimensional surfaces in Mr. Usual excitations 
in wave theory propagate along one-dimensional trajectories. In the 
present paper the role of the medium of propagation of excitations is 
£sj ■ played by the infinite dimensional space of (n — l)-dimensional surfaces 

in Mr. The role of rays is played by n-dimensional solution surfaces 
of a variational problem. Such a generalization of wave theory can be 
\ useful in quantum field theory [2]. 

The paper consists of three sections. In §1 the formula for variation 
of action given by an n-dimensional integral is recalled. In §2 analogs 
of the Hamilton-Jacobi equation and of the canonical Hamilton equa- 
tions are deduced. Besides that, the analogs of the canonical Hamilton 
equations are identified with the equations of characteristics for the 
o analog of the Hamilton-Jacobi equation. The analog of the Hamilton- 

Jacobi equation for n-dimensional variational problems has been found 
in particular cases by many authors, see, for example, [3] and references 
therein. The book [3] also contains a theory of integration of the analog 
of the Hamilton-Jacobi equation, but the equations of characteristics 
are written in [3] in a form somewhat different from that of the present 
paper. In §3 an analog of the wave Schrodinger equation is deduced, 
and an analog of the transport equation is written out in the first order 
of quasiclassical approximation. Besides that, the Cauchy problem for 
the analog of the Hamilton-Jacobi equation is formulated and solved. 

I. The formula for variation of action 

1.1. The main formula. Consider the variational problem for the 
integral 
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Here x 1 , . . . , x n are independent variables, z 1 , . . . , z m are dependent 

7 i dz 1 



variables, z % } = Jpj, and integration goes over an n-dimensional surface 



D with boundary dD in 
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The Euler-Lagrange equations for the functions z l (x) read 



(2) 
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l,...,m. 



Here F z i and F z i denote the partial derivatives of the function F with 
respect to the corresponding arguments. 

Assume that for each (n — l)-dimensional parameterized surface C 
in M m+n given by the equations 
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and sufficiently close to a fixed (n — l)-dimensional surface, there exists 
a unique n-dimensional surface D with boundary 91? = C which is an 
extremal of the integral (1), i.e., the graph of a solution to the Euler- 
Lagrange equations (2). Denote by S = S(C) the value of the integral 
(1) over the surface D. The goal of this section is to deduce the well 
known formula for variation of the functional S. 

To this end, let us represent the integral (1) in the parametric form 
with the parameters s 1 , . . . , s" _1 , t, assuming that the boundary dD 
corresponds to one or several fixed values of the variable t. We obtain 

J= f F{x 1 (s,t),...,x n (s,t),z 1 (s,t),...,z m (s,t), 
Jd 

d(z\x 2 ,...,x n ) d(x\...,x n ) 
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d{x 1 ,...,x n ~ 1 ,z m ) d{x 1 . 
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dfos 1 ,...,*"" 1 ) d(t,s\...,s n - 1 ) 
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ds 1 . . . ds n 1 . Denote 



is the Jacobian, ds 

the integrand in the latter integral by z\ x J gk , z l gk ,xl, z\). Then we 

have (omitting the summation sign over repeating indices) 

(5) 

5J= (<$> x i5x j + QjSz* + <f> x j 5x j gk + <5> x j5x{ + $ 4 5^* + $ 4<^t) dtds 
Jd » k 1 Sfc 



D 



9 * d * 



<9s fe 



at 



+ <E>i - — — r$~i - -=-<&,i 5-2 



at"- 

dtds. 



dtds 



EXCITATIONS PROPAGATING ALONG SURFACES 



3 



It is not difficult to see that the first integral here vanishes due to the 
Euler-Lagrange equations (2). And the second integral equals, by the 
Stokes theorem, 

(6) 5S = J (® x j6x j + <S> z i6z^ ds. 

Hence 



(7) 



where 



d(x 1 , . . . , x l , . . . , x r ' 



Here the cap over a variable means that the variable is omitted; the 
summation sign over the index % repeated twice is omitted. Note that 
in the formula for ffi , the coefficients before the Jacobians coincide 
up to sign with the components of the energy-momentum tensor. Also 
note that the quantities p % and W depend on the numbers x J t , z\ only 
through the numbers z l xj characterizing the tangent space to the n- 
dimensional surface D. 

Thus, we obtain the formula 

(9) 5S= I (J2p i5zi ~J2 HjSxj ) ds 

analogous to the formula for the differential of the action function in 
classical mechanics. This formula easily implies the Noether theorem, 
as in classical mechanics. 

1.2. The properties of the quantities p l and HK It is easy to see 
that the quantities p % and H j satisfy the relations 

(10) p l zl k - Wx j sk =0, fc = l,...,n-l, 
and 

(11) p l z\ - H j 4 = $. 
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The latter equation is just the Euler theorem on homogeneous func- 
tions, since the function $ is homogeneous of degree 1 in the variables 

Hi x t- 

Let us call by the tangent element the data consisting of an (n — 1)- 
dimensional surface C and a family of n-dimensional planes tangent to 
the surface at each point of the surface. These planes can be charac- 
terized by mn variables z l xj satisfying m(n — 1) equations 

(12) 4,-a£ fc = z' fc> fc = l,...,n-l. 

On the other hand, these planes or, more precisely, the dual conormal 
lines can be characterized by m + n variables p\ H j subject to (n — 1) 
equations (10). One passes from the former set of variables to the latter 
one by the duality transform (8). 

A count of parameters shows that the quantities p\ should satisfy 
one more equation. To write it explicitly, assume for a moment that the 
function $ of the variables x\ , z\ is convex, so that it yields a norm on 
the quotient space of the space of vectors (z\, x{) by the transformations 

(13) ^ 

z t ~^ z t + / y a k z l k 

with arbitrary a k . Denote by H the norm on the dual space of vectors 
(p\H j ) satisfying relations (10), i.e., for arbitrary (p\ H j ) satisfying 
(10) put 

(14) H(p\-H j )= max p^-H'^. 

*(4,xJ)=l 

Then it is not difficult to deduce the following relation between the 
quantities p\ H j (8): 

(15) H(p\ -H j ) -1=0. 

The function 7i is determined as a function of m + n variables not 
uniquely, but only up to adding a product of any function with the 
left-hand side of one of equations (10). In the general case of not 
necessarily convex function $, let us mean by equation (15) the n-th 
equation for the quantities p\ HK 

Using the n equations (10) and (15), one can (in general) express the 
quantities H j as functions of p % (and of x l , z\ x l Sk , z l Sk ): 

(16) W = W(x l ,z\x l sk ,zl k ,p l ), j 1 //. 

The change of variables (8), from mn variables z l xj subject to m(n— 1) 
equations (12) and from the function F to m variables p % and to n 
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functions W (16), can be called the modified Legendre transform. Up 
to notations it coincides with the usual Legendre transform. Indeed, if 
instead of the variables z% we consider m variables z\, then the function 
H^xl of the variables p 1 will be the (usual) Legendre transform of the 
function <£>, see formulas (8) and (11). Therefore, since the Legendre 
transform is a duality, we obtain the formulas 

(17) z\ = 4,4 = (tfxbp, = n y,. 

Besides that, if we differentiate equations (10) with respect to the vari- 
able p\ then we obtain the system of equations 

(18) z\ k = zi jX { k = Hy gk , k = 1, . . . , n - 1 
which, together with (17), gives 

(19) 4=#j- 

1.3. Geometric optics in the space of (n — l)-dimensional sur- 
faces. The above formulas become more clear if one has in mind the 
following geometric picture. Consider the space whose points C are 
(n — 1) -dimensional parameterized surfaces (3) in M m+n . Consider 
propagation of light in this space by the Fermat principle, the time 
of propagation of light along a trajectory C(t), to < t < ti, being given 
by the n-dimensional integral (4). Let us call by the light rays the tra- 
jectories C(t) for which integral (4) takes a stationary value. Here one 
has the following essential difference with usual geometric optics: there 
are not one but infinitely many rays passing through two given points. 
These rays correspond to various parameterizations of n-dimensional 
solution surface D of the variational problem. 

Let us fix a point C , and denote by S(C) the time of propagation of 
light from the point C to the point C. Let us call a level hypersurface 
S(C) = T of the function S by a wave front emitted by the point Co 
at the time T, and denote this front by W(T). Then it is natural to 
assume that the following Huyghens principle holds: the wave front 
W(T + AT) is the common tangent hypersurface to the fronts emitted 
by the points of the wave front W(T) at the time AT (which may 
be called spheres of radius AT). Here again one meets the following 
difference with usual geometric optics. At a given point of the front 
W(T + AT), there are infinitely many spheres tangent to the front. 
These spheres correspond to partitions of the surface D of "square" 
T + AT into two pieces of "squares" T and AT. Similarly, each sphere 
is tangent to the front W(T + AT) at infinitely many points. 

Now if we take AT infinitely small, then we obtain the following 
picture. Consider a sphere of an infinitely small radius dT with the 
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center at a point C of the front W(T). Then the function S achieves 
its maximum at the intersection points of the sphere with the rays 
passing through the points Co and C . Hence the tangent hyperplane 
to the wave front at the point C is parallel to the tangent hyperplane 
to the sphere at any of these intersection points. In particular, the 
slope of the tangent hyperplane does not depend on a concrete point 
of the sphere, and depends only on the tangent planes to the surface D 
(the tangent element, see Subsection 1.2). From these considerations 
one can deduce formulas (7,8). 

2. Analogs of Hamilton's formulas 

2.1. Analog of the Hamilton-Jacobi equation. Substituting (7) 
into equations (10,15) or into (16), we obtain 

SS , 5S 



(20) 



\8j(s)' 5xi(s) 1 ' 



or 



(21) ^ +H i^ z '^' zU 'mt) =0 ' J ' " 

The system of equations (20) or (21), relating the variational derivatives 
of the functional S at the same point s, can be naturally called the 
analog of the Hamilton-Jacobi equation. The first n — 1 equations of 
the system (20) correspond to the fact that the function S does not 
depend on the parameterization of the surface C . 

Example (the scalar field in two-dimensional space-time) [1]. Let 
F(x, y, z, z x , z y ) = \{z 2 x - Zy) +p(z), where p(z) = ^z 2 + ... is a poly- 
nomial in z. A calculation gives the following analog of the Hamilton- 
Jacobi equation: 

5S 5S 5S n 
ox dy oz 

(22) ' ' f$S\ 2 ,\ , , .,, , ^ 5S 5S n 

2.2. Analogs of the canonical Hamilton equations. Let us as- 
sume again that the surface D is parameterized by the coordinates 
si, . . . , s„_i, t, and let us find the equations expressing the dependence 
of the variables p\ z % on t, assuming that the dependence of x^ on (s,t) 
is given. 
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To this end, consider the integral 
(23) I [ {p l z\ - H j x 3 t ) dsdt 



as a functional of functions p\ z\ in which IP is the function (16), and 
let us write the Euler-Lagrange equations for this functional. Varying 
p l , we obtain equation (17). It implies that the function p % z\ — Wx\ 
is the Legendre transform of the function Wx\ with respect to the 
variables p % . Hence it coincides with the function $ defined in §1. 
Hence the Euler-Lagrange equations for the integral (23) are equivalent 
to the Euler-Lagrange equations for the function $ and therefore to the 
Euler-Lagrange equations (2). It remains to vary the integral (23) with 
respect to z\ which yields the following form of the field equations: 

5 



(24) " ~W 



H 3 x 3 t {s) ds, 



These equations can be naturally called analogs of the canonical Hamil- 
ton equations. 

2.3. Analogs of the canonical Hamilton equations as the equa- 
tions of characteristics for the analog of the Hamilton-Jacobi 
equation. The analog of the Hamilton-Jacobi equation can be writ- 
ten in the form (9), where are the functions (16). In this form, the 
analog of the Hamilton-Jacobi equation is a system of equations for 
the coordinates p l (s), fP'(s) of the tangent hyperplane to the graph of 
a solution S. At a fixed point (x J '(s), z l (s), S) of the graph, this system 
gives a family of tangent hyperplanes parameterized by functions p l (s). 
The common tangent hypersurface to this family is the analog of the 
Monge cone (see [4], Ch. II). It is given by the equations obtained by 
taking the variational derivative of equation (9) with respect to p l (s), 
i.e., 

(25) zl(s) = Hy t (s), 

where z\ = x 3 t = These are equations (17) forming the first 
half of the analogs of the canonical Hamilton equations. 

Let us call a 4-tuple of functions (x J '(s), z' t (s),p t (s), W(s)), where 
Hi(s) = — g^wpj , p l (s) = fi^fn* ; by an integral element of a solution to 
the analog of the Hamilton-Jacobi equation. An integral element gives 
a point C and the coordinates of the tangent hyperplane to the graph 
of the solution at this point. Assume we are given a one-parameter 
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family of integral elements (x^(s, t), z l (s, t),p l (s, i), H*(s, t)) of a given 
solution S. Let us call the family focal (as in the usual theory of the 
Hamilton-Jacobi equation) if the analog (25) of the Monge equations 
holds. Then one must also have the equations 
(26) 

t _ d SS f u 5 2 S f { , 5 2 S , 

Pt ~ dtdz^s) ~ J XAs) 5z%s)5xi{s') + J Zt{S) 5z*(s)5z*'(s>) 

= - J 4{s')j^ ) HKs'J{z,x))ds' + J H^Wj^/Wds* 

where in the last integral the functions fP are considered as functions 
of independent variables z\p\ The latter computation is somewhat 
formal but it is not difficult to make it rigorous. 

Thus, a focal family of integral elements belongs to the graph of a 
solution S only if equations (25,26), i.e., the analogs of the canonical 
Hamilton equations, hold. If we pass from integral elements (x J , z l , p l , W) 
to tangent elements (x J , z\ z l xj ), then the last condition means that this 
family of tangent elements covers a surface z = z(x) which is the graph 
of a solution to the Euler-Lagrange equations (2). So the analogs of 
the canonical Hamilton equations can be naturally called the equations 
of characteristics for the analog of the Hamilton-Jacobi equation. And 
solution surfaces of the Euler-Lagrange equations (2) can be called 
analogs of rays. 



3. Analog of the Schrodinger equation 

3.1. Analog of the Schrodinger equation. Suppose that the func- 
tions H j (16) are polynomials in the variables p\ Let us make the 
following substitution in the analog (21) of the Hamilton-Jacobi equa- 
tion: 
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Here i is the imaginary unit, h is a very small constant. We obtain a 
system of linear variational differential equations which can be natu- 
rally called the analog of the Schrodinger equation: 
(28) 

- ih sW) + H3 { x ' ^ ^ - ih jJhj) * = °' 3 = 1 -- n - 

Here ty = ty(C) is the unknown complex- valued functional of functions 
xi{s),z l (s); it is assumed that in the polynomials W the operators 
—ih Sz i'r s \ stand to the right of the coefficients of the polynomial. 

The system (28) can be also written in the form of type (20) if we 
assume that the left-hand side of equation (15) is a polynomial with 
respect to the variables H^p 1 : 

Sty { , Sty j _ 
, , S^(s) Z ' k + 6^j X ' k=0, * = 1 --- n - 1 > 

(29) r ( * * \ l 

The first n — 1 equations of the system (29) can be replaced by the 
condition that the functional ty{C) does not depend on the parameter- 
ization of a surface C. 

Example. In the example considered at the end of Subsection 2.1, 
we obtain the following analog of the Schrodinger equation: 
(30) 

1 ( ,2 ^ 2 A , 2 2^ / w , f \ 

2 {- h 6W + **) + ^ - Vl) ~ lK \ Xs W) + J = 

6^ Sty Sty 

s Sx(s) ^ s Sy(s) s Sz(s) 

3.2. The quasiclassical approximation. Let us substitute into the 
system of equations (28) a functional of the form 

, , , , iS(C) 

(31) ty(C) = a(C)e^. 

Equating to zero the principal part with respect to h, we obtain the ana- 
log (21) of the Hamilton- Jacobi equation for the functional S. Equat- 
ing to zero the coefficient before h, we obtain the analog of the transport 
equation: 

< 32 > j^+^y.^y^-o, *=!,..,», 

where p l = s £f^ . This equation can be interpreted in the following way. 
Consider a one-parameter family of integral elements of the solution S 
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of the analog of the Hamilton- Jacobi equation satisfying the equations 
of characteristics, i.e., the analogs of the canonical Hamilton equations. 
Then, as shown above, the corresponding tangent elements cover an 
n-dimensional surface D: z = z(x) which is a solution to the Euler- 
Lagrange equations. Also relation (19) holds. Therefore, the analog of 
the transport equation means that the "amplitude" a(C) is constant 
along the one-parameter family. 

Suppose that the solution to the Cauchy problem for the Euler- 
Lagrange equations (2) with the given tangent element for t = t is 
unique. Then for any other one-parameter family passing through the 
same surface C = C(t ), the corresponding tangent elements cover 
the same surface D. Suppose that the analog of the Monge equations 
(25) has a solution z l (s,t) for any functions x^(s,t) and the initial 
conditions (x^ (s,to), z l (s,to))- Then any (n — l)-dimensional surface C 
on the surface D can be included into a one-parameter family. Hence, 
the integral element at the point C corresponds to the tangent element 
tangent to the surface D, and the value of the function S at the point 
C is given by the formula 



where D\ is the piece of the surface D bounded by the (n— l)-dimensional 
surfaces Cq and C. Formula (33) follows from formulas (9), (11) and 
from the argument of Subsect. 2.2. We obtain that the function a(C) 
does not depend on an (n — l)-dimensional surface C on the surface 
D, but depends only on the surface D itself. This is the analog of the 
statement from quasiclassical approximation of quantum mechanics or 
from short-wave limit of wave optics that the amplitude of the wave 
function is constant along the rays. 

In [3] the sets of integral elements corresponding to all tangent ele- 
ments tangent to the same solution surface D of the Euler-Lagrange 
equations are called the characteristics of the analog of the Hamilton- 
Jacobi equation. Thus, we obtain that in general for a solution S the 
space of surfaces C is covered by a family of pairwise non-intersecting 
projections of characteristics along which the function a(C) is constant. 
The arising family of surfaces D can be called the analog of the notion 
of field of extremals in one-dimensional variational calculus. 

3.3. The Cauchy problem for the analog of the Hamilton— 
Jacobi equation. Let us pose the problem to find a solution S of 
equations (21) which takes the given values S(C ) = U(C ) on the 



(33) 




$ dsdt, 
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surfaces Co of the form 

(34) x j {s) = V{s), z i (s)=a i (s), 

where b*(s) are the given functions, and a l (s) are arbitrary functions. 
To construct solution of this Cauchy problem, note that at each point 
C one can define an integral element p l (s) = ^^y, and therefore a 
tangent element. For a given point C let us pass through this point a 
projection of a characteristic, i. e., let us find Co such that the solution 
surface D of the Euler-Lagrange equations passing through the surfaces 
C and C, would be tangent to the above constructed tangent element 
along the surface C . Assume that such C and D exist and are unique. 
After that let us apply formula (33). 

The thus constructed function S can be obtained also using the con- 
struction of common tangent hypersurface generalizing the Huyghens 
principle (see 1.3). Assume that for each surface C (3) sufficiently 
close to a fixed (n — l)-dimensional surface, and for each surface Co 

(34) sufficiently close to some other (n — l)-dimensional surface, there 
exists a unique n-dimensional solution surface D of the Euler-Lagrange 
equations with the boundary C U C. Denote by J(C , C) the integral 
of the function $ over the surface D. Then for each Co the function 
Sc {C) = I (Co, C) + C/(Co) is a solution of the analog of the Hamilton- 
Jacobi equation (full integral, cf. [1]). The graph of the function S(C) 
is the common tangent hypersurface to the family of the graphs of 
the functions Sc (C) for all Co. Indeed, the common tangent hyper- 
surface is obtained by expressing the quantities Oj(s) as functions of 
xi(s'), z l (s') from the equations 

S(C)=I(C ,C) + U(C ), 

(35) J]_ + JIL = o i = i 

5a l (s) 5a l (s) 

These equations imply that the function S(C) coincides with the one 
constructed above. Besides that, its variation at the point C coincides 
with the variation of the function S Co (C). Therefore, the function S 
satisfies the analog of the Hamilton- Jacobi equation, and its integral 
element at the point C corresponds to the tangent element tangent to 
the surface D. 

The analog of the field of extremals for the solution S is the family 
of solution surfaces D of the Cauchy problem for the Euler-Lagrange 
equations with the given tangent elements along the surfaces Co- The 
characteristics are the varieties along which the graph of the common 
tangent hypersurface touches the graphs of solutions Sc (C) from the 
full integral. 
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Thus, we see that the theory of integration of the Hamilton- Jacobi 
equation from [4], Ch. 2 can be generalized to multidimensional varia- 
tional problems. 
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